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Периодические сигналы широко используются в измерительной технике в качестве тестовых сигналов, в медицине в качестве 
стимулов, в строительной индустрии при акустических испытаниях строительных конструкций, в промышленности при испыта- 
ниях механизмов и машин и т.д. Известны и широко применяются приборы и устройства для их воспроизведения — генераторы 
периодических сигналов. Развиты и методы построения автоколебательных систем генераторов для воспроизведения периоди- 
ческих колебаний простейших форм. Тем не менее, задача синтеза математических моделей автоколебательных систем до кон- 
ца ещё не решена. В статье рассматривается метод стационарных автоколебаний, на основе которого несложно построить сис- 
тему уравнений, решением которой являются колебания заданной формы. 


Введение 

Периодические электрические сигналы с пред- 
писанной (заданной) формой используются в изме- 
рительной технике для питания различных датчиков 
физических величин, измерения передаточных ха- 
рактеристик четырёхполюсников. Они применяют- 
ся в метрологической практике при поверке средств 
измерений параметров переменных напряжений, в 
медицине для кардиостимуляции, в автоматике для 
создания программных движений, в строительстве 
при виброиспытаниях различных конструкций и т.д. 
Поэтому проблема получения сигналов с заданной 
формой или с заданным спектром продолжает оста- 
ваться актуальной. Для решения этой проблемы раз- 
виваются два способа: аналоговый и цифровой. Ана- 
логовый способ основан на формировании сигна- 
лов непрерывной динамической системой, а цифро- 
вой базируется на воспроизведении сигналов диск- 
ретной динамической системой. К достоинствам 


первого способа относят практически неограничен- 
ный диапазон частот, а к преимуществам второго - 
простоту воспроизведения сигналов сложных форм. 
Проблема создания динамической системы с задан- 
ным поведением имеет не только прикладное, но и 
самостоятельное теоретическое значение в теории 
колебаний. Для решения указанной проблемы мож- 
но использовать метод локализации [1], метод уп- 
равления по производной [2] и др. В работе [3] авто- 
ром предложен метод синтеза автоколебательных 
систем с заданной синусоидальной формой колеба- 
ний в стационарном режиме. Этот метод назван ме- 
тодом стационарных автоколебаний. Далее этот ме- 
тод применяется для синтеза систем, воспроизводя- 
щих периодические колебания других форм. 

Постановка задачи исследований 

Метод стационарных автоколебаний основан на 
том, что в динамической системе второго порядка 
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х+Г(х) = $(х,х), 


( 1 ) 


•• </ 2 о 

гдех= ^ 


- линейный дифференциальный опе 

ратор, Р() - нелинейный оператор, # = 

стабилизирующий оператор при определенных ог- 
раничениях на Г и § возможны периодические авто- 
колебания со строго заданной формой. Эти ограни- 
чения сводятся к выполнению следующих равенств 

Хст -г р(х ст ) — 0, ^ 2 ^ 

,?(Іст, Хст) = 0. 


где х ст (/) - заданная форма стационарных колебаний. 

Первое уравнение описывает консервативную 
(по АЛ. Андронову [4]) колебательную систему, а 
второе уравнение - управляющую систему. На рис. 1 
показана блок-схема динамической системы (1). 



Рис. 1. Блок-схема автоколебательной системы 


х{1) 


Верхний блок (консервативная колебательная 
система - ККС) отображает действие первого урав- 
нения системы (2), а нижний (управляющая систе- 
ма -УС) - действие второго уравнения. При таком 
представлении в одной из связей появляется сиг- 
нал который в неявной форме присутствует в ма- 
тематическом описании системы (1). Этот сигнал 
формируется управляющей системой. Если коле- 
бания в консервативной системе находятся вне 
стационарного режима, то управляющая система 
создаёт такое управляющее воздействие и(і) под 
действием которого колебания стремятся к стацио- 
нарному режиму. Как только система входит в ре- 
жим, которому соответствуют периодические ста- 
ционарные колебания х а {1), сигнал управления и(1) 
становится равным нулю, и управление прекраща- 
ется. При таком подходе требуемая форма колеба- 
ний определяется консервативной частью систе- 
мы, а управляющая система обеспечивает лишь ус- 
тановление стационарного режима и парирует воз- 
можные отклонения от него. В этом состоит основ- 
ное и принципиальное отличие предлагаемого ме- 
тода от метода управления по производной и мето- 
да локализации. Ниже проводится результаты фор- 
мализованного синтеза ур. (1). 


Основные результаты 

Рассмотрим возможности реализации колеба- 
тельных систем (1) с предписанной формой коле- 
баний. Первое уравнение системы (2) представляет 


собой уравнение консервативной колебательной 
системы. Для него известно [5], что если существу- 
ет вторая производная х(1) и выполняется условие 
х-Р(х)>0 или х-х<0 при хфО, (3) 

то оно имеет однопараметрическое семейство пе- 
риодических решений х=х„(1). Задача построения 
консервативной части уравнения колебательной 
системы сводится к нахождению выражения для 
функции Р(х). Рассмотрим несколько методов на- 
хождения этой функции. 

В начале рассмотрим метод, основанный на 
представлении Р(х) в виде степенного многочлена 
[6]. Условию (3) удовлетворяют нечётные одноз- 
начные функции Р(х), для которых допустимо раз- 
ложение Дх)=Ёс и _,х“ -1 . При этом периодическое 
решение х„(1) может быть представлено бесконеч- 
ным рядом Фурье. В общем случае полученная 
функция х а (1) может иметь особенности в точках, 
кратных Т/2. Поэтому для построения уравнения 
динамической системы х ст (1) необходимо сглажи- 
вать, например, ограничивая ряд конечным чис- 
лом членов разложения N в ряд Фурье, т.е. 

х\Л0 = ^Ь 2і _ 1 8ІП 

к=1 


2 к 


(2к -1)1 


Тогда, удерживая в выражении для / то же ^членов 
разложения, получим 

,.ч * ^.ч-, ^ , I 2Т1 , 


^(0+^* я (0]=-Е*2*-іЬг( 2 *-і) I х 

к= 1 V 


Х81П 




3 7=1 


уУ- 1 


(0- 


Очевидно, что 

*Г‘(0=«)+*0, 


где К0= Ё 

/(=N+1 

Откуда 


— (2А--1)? 

т 


- остаточный член. 


х ст (0 + Р[х'„ (0] = Ё 5ІП [рг~( 2к ~ !) 1 

к = 1 V ^ 

Х |Х С 27 -іѴі, 27-1 "(у - ) 2 *- 1 )] Ъ 2*-1 і 

+Х<Ѵі Ё ^-,.27-і5ІП^(27:-1)Г 
у=1 *=7Ѵ+ 1 V 


X 


+ 


Если с 2] - 1 определить из условия равенства нулю 
выражения в фигурных скобках, т.е. эти коэффи- 
циенты считать решением системы линейных ал- 
гебраических уравнений следующего вида: 

К- 1 = о, к=йт, (4) 


то получим уравнение 

х ст (о+л* ст (0]=т (5) 

в котором /?( і)=Ъ: 1н {[хТ'(і)-хТ'{і)] тем ближе к нулю, 
чем больше А*, при котором правой частью уравне- 
ния (5) можно пренебречь. Число удерживаемых 
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членов _/Ѵ* выбирается из условия минимизации не- 
вязки между решением однородного уравнения (1) 
при начальных условиях 

х ст (0)=0,х сі (0) = |> 24 

к = 1 

**ст(0=ХХ-і'3 Іп 

к = 1 

При этом Р(х)=^с яа х 21 ‘-' и уравнение консерва- 
тивной части принимает вид 

Ха{1) + '^ І с 1к -і-[х СІ (!)] 2к - 1 20 . ( 6 ) 

к = 1 

В ряде случаев, при сложной форме заданных к 
синтезу колебаний, когда И’ велико, функция Дх) 
может описываться степенным полиномом высокой 
степени, тогда в качестве стационарного колебания 
можно принять первую или вторую производную 
некоторой вспомогательной функции ^.Д), такой 
что требуемая форма определяется либо её первой 

Х а {І)=ІМ), ли бО ВТОРОЙ ПРОИЗВОДНОЙ Х СІ (?) = І'ст(0- 

Рассмотрим метод построения функции Дх) в 
виде разложения по многочленам Чебышева. 

Пусть заданное колебание содержит только чёт- 
ные гармонические составляющие, т.е. является 
решением вида 

N ~ 

х„(!) = ^а к -со8 {Щ-кі). 

к = 1 1 


2 л 


(2А--1) 


2л 


(2к-\)1 


Подставляя это решение в первое уравнение 
системы (2), получим 


-^(^к) 2 а к -со&(^кІ) + Р 



= 0 . 


Учитывая, что тригонометрические функции 
можно представить через многочлены Чебышева, 
имеем уравнение, 

ю 2 к 2 а к -Т к (х) + Р ^ а к -Т к (х ) 

к = 1 |_&=1 


= -п(х)+Р[ѵ(х)} = 0, (?) 


2л 

где х = С05—1 = С08 Ш. 

Полученное уравнение по известным Дх) и ѵ(х) 
позволяет определить функцию Р, например, пос- 
редством интерполяции степенными многочлена- 
ми или сплайнами зависимости ѵ/{х)=Р[ѵ(х)\. 

Теперь рассмотрим второе уравнение системы 
(2), которое представляет собой уравнение управ- 
ляющей системы. Найдём его вид в общем случае 
из уравнения консервативной части, для чего 
представим его в виде системы. 

\У = -Р(х), 

\х = у. 

Как известно [4], эта система имеет первый ин- 
теграл 

С = ^ + С(х), (8) 


где С(х) = С - постоянная интегрирования. 

0 

Известно также, что первый интеграл остаётся 
постоянным вдоль каждого из решений системы. 
На фазовой плоскости первому интегралу уравне- 
ния консервативной системы соответствует семей- 
ство (континуум) вложенных друг в друга замкну- 
тых линий (траекторий). Для этого функция С{х) 
должна иметь глобальный минимум [4] . При этом 
замкнутым линиям соответствует особая точка ти- 
па "центр". В континууме траекторий одна из них 
соответствует заданному стационарному периоди- 
ческому колебанию х с Д)- Отклонение от этой тра- 
ектории (при определённом С) 1 или иначе переход 
на другую траекторию меняет значение первого 
интеграла. Фактически константа С задаёт пара- 
метры колебания: его амплитуду и частоту. При 
конкретном значении С переход на другую траек- 
торию приводит к нарушению равенства (8). Воз- 
никает разность между левой и правой частью (8). 
Знак и значение этой разности зависит от отклоне- 
ния от заданной траектории, т. е. от стационарного 
колебания х ст (/). Это свойство первого интеграла 
используется для управления колебаниями. Ведь 
ясно, что параметры соответствующих колебаний 
консервативной системы зависят от начальных ус- 
ловий и при малейшем их изменении существенно 
меняются. Для независимости параметров колеба- 
ний от начальных условий или воздействия возму- 
щений и превращения их в автоколебания необхо- 
димо особую точку системы "центр" превратить в 
"фокус", причём так, чтобы на стационарной траек- 
тории ур. (8) превращалось в тождество. Для прев- 
ращения "центра" в "фокус", введём инкрементную 
составляющую. Таким образом, дифференциаль- 
ное уравнение метода стационарных автоколеба- 
ний можно представить в виде 

х + Р(х) = е х[С-(х) 2 /2-С(х)], (9) 

где е - коэффициент (малый параметр). 

Это уравнение позволяет строить динамические 
автоколебательные системы при заданной форме ко- 
лебаний х ст (/), заданной аналитической функции Р(х) 
или при заданном спектральном составе колебания. 

Практика применения изложенной методики к 
синтезу колебательных систем ( 1 ) показывает, что усло- 
вия (3), налагаемые на функцию Дх), являются очень 
сильными и приводят к существенным ограничениям 
класса воспроизводимых колебаний. Например, если 
заданное к синтезу колебание имеет на половине пери- 
ода значения второй производной равные нулю при 
хфО, то такие колебания не удовлетворяют ур. (3). Нап- 
ример, при построении колебательной системы для ге- 
нерирования колебаний, содержащих первую и третью 
гармоники х ст (?)=а $іп(0+Д 8т(ЗД, воспроизводимыми 
могут бытъ колебания удовлетворяющие неравенству 
Р/а< 0,3. По существу, неравенства (3) ограничивают 
класс воспроизводимых колебаний такими, которые 
имеют на фазовой плоскости только один центр. Одна- 


‘Константа С находится путём интегрирования первого слагаемого (8) за период колебаний. 
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ко изложенную методику можно распространить и на 
случай, когда особых точек типа центр может быть 
больше одной и, соответственно, класс воспроизводи- 
мых колебаний может быть существенно расширен. 
Важно только, чтобы на самом предельном цикле не 
было особых точек и, чтобы стационарная фазовая тра- 
ектория охватывала бы все частные замкнутые циклы, 
соответствующие особым точкам. 

Ниже приведены нетривиальные примеры при- 
менения уравнения (9) для синтеза автоколеба- 
тельных систем, подтверждающие сказанное. 

Пусть к синтезу задано колебание, х„(і)= 
=-1,4сс 8Іп(/)+3,2 р 8Іп 3 (7) график которого показан 
на рис. 2, а в виде линии х(1), а его вторая производ- 
ная - на рис. 2, 6. Как видно, вторая производная 
этого колебания равна нулю при хфО в точках і=п/ 3, 
/=2я/3 и др. При этом не выполняется условие (3). 
Более того, функция Р(х) оказывается неоднознач- 
ной. Общие вопросы построения автоколебатель- 
ных систем с такими функциями ещё не решены. 

Поэтому примем, что х а (1)='і{і), где ф) - неко- 
торое вспомогательное колебание (на рис. 2, а оно 
показано пунктиром). Тогда, выполняя все опера- 
ции синтеза, получим по формулам (4) выражение 
для функции Т(г)=-1,9г+9,57г 3 +21,32^ 5 +44,02г 7 - 
-64, 05^+57, 89г п -28,76г 13 +5,96г 15 . 



- 10 

0 і 2-я б ) 

Рис. 2. Г оафики а) заданного к синтезу сигнала и б) его вто- 
рой производной 

На рис. 3 показан график функции Р(т) (рис. 3, а), 
а так же график погрешности воспроизведения коле- 
баний (рис. 3,6). Видно, что погрешность не превыша- 
ет 1 %, что считается приемлемым для задач синтеза. 

Фазовые портреты возбуждения автоколебаний 
представлены на рис. 4. 

Видно, что в данной системе имеет место слу- 
чай двух особых симметричных точек типа "центр". 
Видно так же, что введение в управляющее воздей- 
ствие инкремирующе - демпфирующих составля- 
ющих привело к превращению "центров" в "фоку- 


сы", в результате чего колебания из областей осо- 
бых точек приходят на стационарную траекторию. 



-2 г 2 а ) 



Рис. 3. Г оафики: а) функции Р(і) и б) погрешности воспро- 
изведения колебаний 



-1.5 7.(1) 1.5 а ) 



- 1.5 7 .( 1 ) 1.5 б) 

Рис. 4. Фазовые портреты установления автоколебаний при 
начальных условиях а) внутри и б) вне стационарно- 
го предельного цикла 

Приведём ещё один пример применения метода 
стационарных автоколебаний к синтезу с заданной 
функцией Р(х), например, в виде 

Р(х) = 16х 5 - 20х 3 +5х= Т 5 (х), 

описываемой многочленом Чебышева пятого по- 
рядка. 
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В этом случае ур. (9) принимает следующий вид 
х + 1 6х 5 - 20х 3 + 5х = 


= е-х\ 


1-(х) 2 



( 10 ) 


1.5 

х(1) 

у(0 

-1.5 

О I 20 а ) 

1.5 

у4(0 
уЗ (О 
у2(0 
уі(0 


- 1.5 

-1.5 х4{1) ,хЗ(1) ,х2(1) , хі(і) 1.5 б) 

Рис. 5. Г оафики: а) стационарного решения уравнения ( 10) и 
6) фазовые портреты при разных начальных условиях 

На рис. 5 представлено решение ур. ( 10) и фазо- 
вые траектории при разных начальных условиях. 
Из графиков видно, что здесь имеет место случай 
трёх "центров'. Тем не менее, все траектории после 
переходного процесса выходят на стационарную 
траекторию - предельный цикл. В зависимости от 
постоянной С и начальных условий здесь возмож- 
ны траектории, охватывающие каждую из особых 
точек или все три. Уравнение метода стационарных 
автоколебаний (9) позволяет строить динамичес- 
кие системы и с заданным спектральным составом. 

Приведём пример синтеза автоколебательной 
системы для воспроизведения колебания с задан- 
ным спектром, например, содержащего первую, 
вторую и третью гармоники с амплитудами а к рав- 
ными единице 

3 

Хст(() = ' со $(кг) = і(Т). 

к=\ 





2 2 а ) 



Рис. 6. Г оафики: а) заданного к синтезу сигнала содержаще- 
го три первые гармоники с амплитудами а 1 =а 2 =з 3 =7 и 
б) функции Р 
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Рис. 7. Функция потенциальной энергии а) и фазовый порт- 


Для такого колебания функции м/(х) и ѵ(х) име- 
ют вид 

ѵг(х) = й 3 (4х 3 -Зх) +а 2 (2х 2 -1) +а г х = %({), 

4 ѵ 3 — 9 г 2 — 1 

ѵ(х) = а 3 - + а 2 ^ +аіХ = і{1), где х = сок/. 

На рис. 6 приведены графики заданного к син- 
тезу сигнала и функции Дх). 


На рис. 7 показан график функции потенциаль- 
ной энергии С(7.) и фазовый портрет системы. Из 
графика видно, что функция (7(Д имеет глобаль- 
ный минимум и, следовательно, консервативная 
система должна иметь континуум циклов. Фазовый 
портрет позволяет сделать вывод о том, что в дина- 
мической системе (9) с нелинейной функцией ЯД, 
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показанной на рис. 7, б, один из циклов превраща- 
ется в предельный цикл системы и все траектории 
стремятся к нему. 



Рис. 8. Спектр колебаний в системе 

На рис. 8 показан спектр автоколебания. Видно, 
что спектр, также как и заданный к синтезу, содер- 
жит три первые гармоники, амплитуды которых 
практически равны единице с погрешностью, не 
превышающей 4 %. Эта погрешность, а также нали- 
чие высших гармоник малой амплитуды при к > 3 
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обусловлены неточной аппроксимацией функции К 

В частном случае, когда заданная форма сигна- 

2п 

ла синусоидальная, т.е. х ст (г) = ып — / = ып ы о і, 

функция Р(х)=а>1х, тогда ур. (9) превращается в 
уравнение 

х + со„х = е й) 0 і[1 - ( со 0 х ) 2 - х 2 ], 

которое имеет решения в виде синусоидальных ко- 
лебаний. 

Заключение 

Таким образом, создан эффективный метод 
синтеза автоколебательной системы с заданной 
формой колебаний и заданным спектром. Метод 
нашёл применение при синтезе автоколебательных 
систем генераторов электрических сигналов. Его 
можно применить при синтезе стимулов в медици- 
не, программных движений механических систем и 
в других случаях. 
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